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■ Abstract. We review the problem of heavy-quark diffusion in the Quark-Gluon 

ON \ Plasma and its ramifications for heavy-quark spectra in heavy-ion collisions at RHIC. 

In particular, we attempt to reconcile underlying mechanisms of several seemingly 
different approaches that have been put forward to explain the large suppression and 
elliptic flow of non-photonic electron spectra. We also emphasize the importance of 
a quantitative description of the bulk medium evolution to extract reliable values for 
Oh, the heavy-quark diffusion coefficient. 
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>- ■ 1. Introduction 

o 

\ The heavy-quark (HQ) mass, toq ^> T c , brings a large scale into the problem of probing 
strongly interacting matter in the vicinity of its (pseudo-) critical temperature, T c . 
This implies that in high-energy heavy-ion collisions charm and bottom quarks (Q=c, 
b) are mostly produced upon initial impact of the incoming nucleons, on a timescale 
r prod ~ 1/rriQ < 0.1 fm/c, while their abundance is expected to be frozen thereafter 
(which is supported by data PQ). Therefore, changes in HQ momentum spectra due to 



O 



reinteractions in the produced medium become an excellent (since rather direct) probe 
of the latter. The large HQ mass also implies that (the approach to) thermalization is 
delayed thus providing a more sensitive means to study the in-medium interactions 
responsible for equilibration. The (approximate) "factorization" of QQ production 
and rescattering renders HQ observables a well-defined probe over the entire range of 
transverse momentum (pt), enabling a comprehensive study of (a) thermalization at 
low pt, (b) a kinetic regime at intermediate pt and (c) "jet-quenching" at high p?. The 
conservation of individual charm and bottom quantum numbers allows to further test 
quark coalescence mechanisms in the hadronization transition. Theoretically, tuq ^> T 
opens the possibility to describe (low-pr) HQ motion in the QGP within a diffusion 
equation [2] which facilitates the extraction of pertinent transport coefficients from 
heavy-ion data. The transition from the elastic to the radiative scattering regime 
is, of course, a key issue which is closely related to items (a)-(c) above. Finally, the 
availability of rather accurate lattice QCD (1QCD) computations of HQ free energies at 
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finite temperature has spurred the hope of being able to define and extract interaction 
potentials which may not only govern the properties of heavy-quark (QQ) bound states 
but also heavy-quark diffusion. The nonperturbative nature of the latter has been 
exhibited in Refs. [21 H] (see Ref. [5] for a recent review). 

Our article is organized into a discussion of HQ diffusion coefficients (Sec. EJ, their 
applications to HQ spectra and observables at RHIC (Sec. [2]) and conclusions (Sec.Hj). 



2. Heavy Quark Diffusion in the Quark-Gluon Plasma 

Brownian motion of a heavy quark in a fluid of light partons is described by a Fokker- 
Planck (FP) equation, schematically written as 

dt dp dp 2 ' 

which follows from the Boltzmann equation to second order in the momentum transfer, 
k, to the heavy quark. The friction (7) and momentum diffusion (D) coefficients, 

1V = J d 3 kw(k,p)k , D = - J d 3 kw{k,p)k 2 , (2) 

are given in terms of transition rates, w, which are typically computed from the elastic 
scattering amplitude of the heavy quark with the (light) medium constituents. Both 
terms in Eq. are essential for the HQ distribution function, fg, to approach thermal 
equilibrium. This feature is highlighted by the Einstein relation, T = D/ (7771,3), which 
also shows that the FP equation involves only one independent transport coefficient. 
The latter is often quoted in terms of the spatial diffusion constant, D s = T/^juiq). 

Early estimates of HQ diffusion in the QGP have utilized "leading-order" (LO) 
perturbative QCD (pQCD) which is dominated by t-channel gluon exchange with 
thermal partons (quarks and especially gluons) [21 El El E] . The infrared divergence for 
forward scattering is regulated by a thermal Debye mass, ~ gT, resulting in a total 
cross section o~q p ~ a s /fi 2 D . However, due to the forward-angle dominated scattering, 
the transport cross section is much smaller; at temperatures around T=300 MeV and for 
a strong coupling constant, a s =0A, the thermal relaxation time for c-quarks amounts 
to r = 7 _1 ~15-20fm/c (larger by a factor of ~mf,/m c for b quarks). This is much 
longer than the expected QGP lifetime at RHIC and would therefore lead to small 
modifications of HQ spectra. 

In Ref. [6] it was suggested that mesonic resonances in the QGP could be operative 
in thermalizing heavy quarks. Pertinent effective Lagrangians where constructed 
utilizing HQ effective theory with a chirally symmetric set of D- and B-\ike resonances 
for c-q and b-q scattering. Within an estimated range of the underlying 2 parameters 
(coupling constant and resonance mass), the HQ relaxation times where found to 
be reduced by a factor ~3 compared to the LO-pQCD results. Pertinent Langevin 
simulations within an expanding QGP at RHIC, augmented with quark coalescence at 
T c , resulted in a reasonable description [9] of non-photonic electron (e 1 * 1 ) observables at 
RHIC [H [in] (cf. left panel of Fig. [5] below). The assumption of resonances calls for a 
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Figure 1. Heavy-quark internal energies as extracted from fits [16j [13 to quenched 
(left) and Nf—2 (right) free energies computed in thermal lattice QCD [18] . 

microscopic treatment to eliminate the free parameters of the effective Lagrangian and 
to assess their temperature dependence. In Ref. [3], a T-matrix equation for heavy-light 
quark scattering was set up, which, after partial wave expansion (L=0,1,. . . ), takes the 
form 

Tl(E cm] q, q') = V£(q, q') +Jdkk 2 V£(q, k) G Qq (E cm , k) T a L {E cm - k, q') ; (3) 

Gq q denotes the intermediate 2-particle propagator which in principle contains in- 
medium selfenergies (mass corrections and width) of the individual heavy and light 
quark. The main idea is to implement the interaction kernel, V£, in a parameter- free 
way by utilizing first-principle information from 1QCD computations in terms of the HQ 
free energy at finite temperature, 

F QQ (T,r) = Uqq(T, r) — TSqq(T, r) . (4) 

The same approach has been applied in recent years to compute heavy quarkonium 
spectral functions in the QGP, with fair success in reproducing independent 
1QCD results for Euclidean-time correlation functions [Til [T2| [T3| HI) . In the 
vacuum, Fqq{t) = Uqq{t) closely resembles the phenomenological Cornell potential 
( Coulomb + confinement) which has been very successful in quarkonium spectroscopy 
and is now understood as the low-momentum limit of QCD with heavy quarks (so- 
called "potential QCD") [15]. The appropriate potential definition at finite temperature 
is currently an open question. An upper limit may be obtained by using the (subtracted) 
internal energy, VqQ(r; T)=Uqq(t] T) — f/gg(r=oo; T), which in the quarkonium sector 
provides the largest binding [13]. The subtraction is required to ensure convergence of 
the Fourier transform of the potential in momentum space in Eq. (J3j). Two examples of 
pertinent potential extractions are shown in Fig. [TJ Their application in the T-matrix 
equation ([3]) additionally includes a relativistic correction to simulate color-magnetic 
interactions [TYJ. The Born approximation to the T-matrix, T£ = V£, recovers the 
results from LO-pQCD within ~10% above T cm =4 GeV (for the same value of a s ). The 
different color channels, a, are accounted for by Casimir scaling of the potentials. It 
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Figure 2. Charm-quark friction coefficients 7 in the QGP. Left panel: 3-momentum 
dependence at 3 temperatures (color code) for: LO-pQCD with fixed a s —0A and 
(j,D=gT (lower 3 lines), heavy-light quark T-matrix plus LO-pQCD for gluons (middle 
3 lines) [4], and pQCD with running a s and reduced infrared regulator (upper 3 
lines) [19]. Right panel: temperature dependence of 7 for LO-pQCD, T-matrix plus 
LO-pQCD (gluons only), pQCD with running a s , and from AdS/CFT correspondence 
matched to QCD [22] with (7=1.5-2.6 [23] in Eq. ©. 



turns out that the resulting in-medium T-matrices are dominated by 5-wave scattering 
in the attractive color-singlet (meson) and -antitriplet (diquark) channels, supporting 
resonance-like structures close the qQ threshold up to T~ 1.7 T c and 1.4 T c , respectively. 
Close to T c , the corresponding friction coefficients are a factor ~3 larger than LO-pQCD 
(see Fig. [2]), but they decrease as temperature increases due to the weakening (color- 
screening) of the lQCD-based potential in the left panel of Fig. [1] (for the potential in 
the right panel of Fig. [1] a slight increase of 7(T) is found). It is interesting to note 
that the collisional dissociation of D and B mesons in the QGP evaluated in Ref. [20] is 
based on a similar Cornell- type potential as the T-matrix [4], but without the inclusion 
of medium effects in the potential. The effect of a reduced formation time of B relative 
to D mesons [20] is encoded in the T-matrix calculations as a mass effect leading to a 
stronger binding (analogous to the heavy quarkonium sector [14J). 

Perturbative evaluations of HQ transport have recently been revisited in Ref. |19j . 
where it is argued that the infrared regulator in the t-channel gluon propagator becomes 
operative at a significantly softer scale than the Debye mass. In addition, a running of 
the strong coupling constant to small scales is implemented. The combined effect on the 
HQ friction coefficient is an approximately ten-fold increase over the Born result with 
Debye-mass regulator, cf. Fig. [2] Under these circumstances the perturbative treatment 
should be revisited and presumably augmented by resummations. 

HQ diffusion has also been evaluated in conformal gauge theories where 
nonperturbative results can be inferred using a correspondence to string theory 
( AdS / CFT) [21] . Defining the friction coefficient as the (inverse) timescale of momentum 
degradation, dp/dt = —•yp, cf. Eq.([T]), one finds for a Af=4 Super- Yang-Mills plasma: 

tt\/AT| ym T 2 
7AdS/CFT = — <-> 7qcd = k • (oj 
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Figure 3. Spatial diffusion coefficient, D s —T/(jmQ), for c (left) and b quarks (right) 
in a QGP for: LO-pQCD with fixed a s =0.4 (dashed lines), effective resonance model 
+ LO-pQCD (bands for T D>B =0.4-0.75 GeV) [6j, T-matrix approach + LO-pQCD 
(gluons only) [A\ and pQCD with running a s (dash-dotted line) [15]. The AdS/CFT 
result, Eq. ((5]), corresponds to 2irTD s =2ir/C~l .5-4 (not shown in the plots). 



The relation to a QCD plasma (second expression in Eq. requires a careful trans- 
lation of the temperature scale and coupling constant; matching the energy densities of 
the SYM plasma and QGP to identify T and using the static HQ force from 1QCD to 
identify the coupling strength, A, one finds [22] C~1.5-4. Note that this identification 
makes direct contact with the 1QCD HQ free energy, as in the T-matrix approach. 

The right panel of Fig. [2] compiles the temperature dependence of the friction 
coefficient, 7(p=0), in the approaches discussed above. The pQCD+running-a s and 
AdS/CFT results are quite comparable but significantly (much) larger than the ones 
from the T-matrix (LO-pQCD), except close to T c . However, it might be in the vicinity 
of T c where the QGP is most strongly coupled providing favorable conditions for the 
validity of the strong-coupling limit. The T-matrix approach exhibits the weakest 
temperature dependence of the displayed curves, with a moderate uncertainty induced 
by different version of the internal energy (however, using the free energy as potential 
would give results closer to LO-pQCD). 

Finally, we compare in Fig.[3]the spatial diffusion coefficients (in units of the thermal 
Compton wavelength, l/(27rT)) at zero 3-momentum (which may be thought of as being 
proportional to the ratio of shear viscosity to entropy density) for charm and bottom 
quarks. All approaches give results fairly independent of temperature and HQ mass, 
except for the lQCD-based T-matrix calculation where the increase with temperature 
indicates a significant loss of coupling strength in the QGP. 

3. Heavy-Flavor Spectra at RHIC 

Several groups have applied the Brownian motion framework to simulate HQ diffusion 
in Au-Au collisions at RHIC, using hydrodynamic [3, [T9J, [23] or expanding fireball [9] 
models. Obviously, a realistic description of the bulk medium (temperature and flow 
evolution) is mandatory for a quantitative extraction of diffusion coefficients from data. 
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Figure 4. Time evolution of nuclear modification factor (central collisions, left panel) 
and elliptic flow (semicentral collisions, right panel) in a QGP fireball at RHIC. [5]. 



Let us first address the time evolution of HQ p r -spectra, characterized by the 
nuclear modification factor Raa(pt) — {.oINaa/ 'dpr) 1 1 (N co udN pp / 'dpx), and elliptic flow 
coefficient, i^Cpt)- The results in Fig. HJ computed in a thermal fireball evolution with 
HQ transport based on resonance+LO-pQCD interactions [9], suggest that the high-p T 
suppression is built up significantly earlier than the elliptic flow. The former feature is 
quite reminiscent to a recent analysis [24] of the empirical system-size dependence of 
high-pr hadron suppression, arguing that parton energy-loss is predominantly operative 
in the first 2-3fm/c. On the contrary, the bulk v-i in hydrodynamic models (used to 
construct the fireball evolution) requires a duration of at least Ar~2-3fm/c to build up 
most of its strength, thus "delaying" the transfer to heavy quarks. 

Next we compare Langevin simulations for charm quarks in semicentral Au-Au 
collisions at RHIC. Representative values obtained for v 2 and Raa in 3 different evolution 
codes (2 hydro, 1 fireball), implementing 3 of the approaches discussed above for HQ 
transport coefficients, are summarized in Tab. HJ Given the complete independence of 



Model [Ref.] 


D s (2nT) 


b [fm] 


,,max 
u 2 


Raa(pt=5 GeV) 


hydro + LO-pQCD [7] 


24 


6.5 


1.5% 


0.7 


hydro + LO-pQCD [7] 


6 


6.5 


5% 


0.25 


fireball + LO-pQCD [9] 


-30 


7 


2% 


0.65 


fireball + reso+LO-pQCD [9] 


~6 


7 


6% 


0.3 


hydro + Eq. © [23] 


21 


7.1 


1.5-2% 


-0.7 


hydro + Eq. © [23] 


2tt 


7.1 


4% 


~0.3 



Table 1. Overview of model approaches (f . column) and input parameters (2. column: 
spatial charm-quark diffusion coefficient, 3. column: nuclear impact parameter) for 
Langevin simulations of charm-quark spectra in Au-Au collisions at RHIC; selected 
values for the resulting elliptic flow (v™ ax ~u 2 (pT=5 GeV)) and nuclear modification 
factor are quoted in columns 4 and 5. 
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Figure 5. Nuclear modification factor (upper panels, central Au-Au) and elliptic flow 
(lower panels, minimum-bias Au-Au) of non-photonic electrons at RHIC [Tl [TU1 125]. 
compared to theory. Left panel pQ: Langevin simulations based on (i) hydrodynamic 
evolution with (upscaled) LO-pQCD HQ interactions (dotted and dash-dotted lines) [7] 
and (ii) expanding fireball with effective HQ resonance interactions (band) [9|, and 
pQCD radiative energy-loss calculations with (upscaled) transport coefficient in a 
static medium [26) : right panel: QGP fireball with T-matrix interactions [4] using 
2 extractions of the HQ internal energy from lattice QCD as potential (see Fig. Q]) 
where the dashed lines do not include heavy- light quark coalescence at T c . 



the calculations, the agreement on the ~30 % level is encouraging (it might improve 
when accounting for finer details, e.g., the results of Ref. [23J correspond to terminating 
the evolution in the middle of the mixed phase while in Refs. [9] it is run to the end 
of the mixed phase, cf. right panel of Fig. H]). 

Finally, theoretical predictions for non-photonic electron spectra are compared to 
RHIC data in Fig. [5j This requires the conversion of the modified c- and 6-quark spectra 
into meson spectra followed by their semileptonic decay. For the radiative energy- 
loss calculations of Ref. [26] and the Langevin simulations of Ref. [7] hadronization 
is treated via independent fragmentation, while in Refs. H] coalescence processes 
with light quarks from the medium are accounted for. The latter enhance the hadron 
(and electron) t> 2 while reducing the suppression (see dashed lines in the right panel 
of Fig. [5]). While the radiative energy-loss calculations (with upscaled pQCD transport 
coefficient) roughly account for the observed suppression, the corresponding elliptic flow 
is underpredicted. This underlines the importance of accounting for the collectivity of 
the expanding QGP medium, transferred to heavy quarks via the diffusion term in 
Eq. ([1]). The current data-theory comparison points at HQ diffusion coefficients in the 
range of D s (2ttT)^A-Q. 
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4. Conclusions 



Several calculations of HQ diffusion in a QGP at moderate temperatures have been 
conducted recently but do not (yet?) show satisfactory agreement, neither in magnitude 
nor in their temperature dependence. We have argued, however, that there is significant 
overlap in the underlying physics, which in most cases is based on a (in-medium) 
color-Coulomb type interaction as implicit in both one-gluon exchange and potential 
models (or in matching to conformal gauge theories). The challenge is thus to reconcile 
the different approaches (perturbative, T-matrix, AdS/CFT, etc.) by revisiting their 
regimes of applicability. The inclusion of radiative processes in a diffusion framework 
remains another challenge. On the empirical front, electron spectra at RHIC indicate 
a small HQ diffusion constant, but also here further scrutiny in the implementation of 
different medium evolution models (including the hadronic phase) needs to be exerted. 
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